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I address here the question of the mutual interplay of strong correlations and 
disorder in the system. I consider random version of the Hubbard model. 
Diagonal randomness is introduced via random on-site energies and treated 
by the coherent potential approximation. Strong, short ranged, electron - 
electron interactions are described by the slave boson technique and found to 
induce additional disorder in the system. As an example I calculate the den- 
sity of states of the random interacting binary alloy and compare it with that 
for non-interacting system. PACS numbers: 71.10.Fd, 71.23.-k, 7 1.55. Ah. 

1. INTRODUCTION 

The list of materials in which carriers are strongly interacting with each 
other and scatter by random impurities is quite long. It comprises inter 
alia high temperature superconducting oxidesEI and various heavy fermion 
alloys.cl Two models are usually used to study such system. It is either 
Hubbard or Anderson model suitably extended to allow for the description 
of real materials. Here the single band Hubbard model has been used, as 
it is the simplest model of correlated system. The disorder is introduced 
into the model by allowing for fluctuations of the local on-site energies £j. 
The main purpose of this work is to study the interplay of disorder and 
correlations in the Hubbard model. I shall present analytical and numerical 
calculations indicating that correlations induce additional disorder in the 
system. Mutatis mutandis disorder in the system affects the parameters 
of the Mott-Hubbard metal-insulator transition. In the recent studies of 
weakly interacting disordered systems it has been found that due to disorder 
the interactions scale to strong coupling limit J! Here to treat many particle 
aspect of the problem we use the slave boson technique, which is known to be 
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qualitatively valid for all strength of interaction. In particular this method 
reproduces the results of the Gutzwiller approximation at the saddle-point 
level. □ 



2. THE MODEL AND APPROACH 

We start here with random version of the Hubbard model 

ija itr i 

The meaning of symbols is standard. The first term describes the hopping 
of carriers through the crystal specified by lattice sites i, j. e« denotes the 
fluctuating site energy - it introduces disorder into the model. U is the 
repulsion between two opnosite spin electrons occupying the same site. In 
the slave boson technique^ one introduces 4 auxiliary boson fields: e^, §i a , 
di such that efii, sf a Si a , dfdi project onto the empty, singly occupied by a 
and doubly occupied site i. The Hamiltonian (1) is in this enlarged Hilbert 
space written as 

H = J2 Uj^JjaC^Cja + Uj^ dfdi + ^(e; - n)4a C i<T , (2) 
ijcr i icr 

where z ia = (1 -dfdi ~sf Sia)~ 1/2 {ef s io - + sj_ a di) ■ (1 -efii - sf- a Si-<r)~ 1/2 ■ 
Equation (2) is strictly equivalent to (1), when the constraints 

ef^i + sfsio- + dfdi = 1 

C fa C ™ = vfaVia + dfdi] <? = |, j (3) 

are fulfilled. In disordered system and at T = OK, the mean field approach 
to slave bosons can be formulated by suitably generalizing the clean limit. 
One introduces the constraints into Hamiltonian with help of Langrange 
multipliers A^ and A^ and replaces boson operators e^, §i a , di by classical, 
site dependent, amplitudes ej, Si a , di. These amplitudes are calculated from 
the, configuration dependent, ground state energy Eqs = (H) by minimizing 
Eqs with respect to all seven parameters e^, Sj CT , di, As a result 



3 2 + Y i _ I /,-.+ ^ _ ,.2 , ,/2 

additional equations which read 



one gets three constraints: ef + s ia + dj = 1, {cf a Ci a ) = s ia + df and four 



A^e; = -~|^Ite \ Y^t ij £ i {cf (T c j(T )t i 
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(u + A« -E>&) * = ~lM f Re E^fe>)^j (4) 

where & = (1 - d? - s^)" 1 /*^ + s ,_ ff d 4 )(l - e? - s^)" 172 - 

The important next step is connected with calculation of (cf a Ci a ) and 
-Eio- = J2j tij£,i( c tr c j<r)£i m * ne presence of disorder. In principle these quan- 
tities do depend on the particular distribution of all impurities (i.e. distri- 
bution of site energies £i). We shall calculate them from the corresponding 
Green's function. We use a versionQ of the coherent potential approximation 
(CPA) to calculate (averaged and conditionally averaged) Green's functions. 
To this end we rewrite the mean field Hamiltonian in the form 



H = Y li jZi\j r i« r i» + J2( £ i - A* + x il ] ) c t c ^ + const • ( 5 ) 

ijcr 



Note, that Hamiltonian (5) contains not only the site dependent parameters 

(2) 

£j but due to correlations there appear also \\J which are the additional 
source of, diagonal and spin dependent, disorder. Besides that, the param- 
eters £j, £j which vary from site to site make effective hopings iy = 
random quantities. This off-diagonal (in Wannier space) disorder is particu- 
larly important. Fortunately multiplicative dependence of effective hopping 
on the kind of atoms at sites i and j is easy to handle within CPA. 

The procedure is standards and one gets the following equation for the 
density of states in paramagnetic phase 

D(s) = -- Im( = C 2 FP(£ + )] y (6) 

vr ! _ [ S(e+) _ {£ _ £t + fl _ A f) )/e 2] F[m] 

and the coherent potential S(e) is determined as a solution of equations 



1 

F[E(z)] = i=Y ^ = ttYGAz) (7) 

Here (• ■ -)i mp means averaging over disorder and e(k) = J2ij Uj e~ lk ( Ri ~~ R ^ 
is electron spectrum in host (clean) material. 
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To solve equations (4) one still has to calculate the quantities Ei a = 
J2j tij£,i{ c ta c ja)ij , which in general depend on the configuration of all impu- 
rities in the system. They can be calculated from the knowledge of the CPA 
Green's functions in the following way. We assume that site i is described 

(2) 

by actual parameters i.e. £», X ia , £j, while all other sites in the system are 
replaced by effective ones described by the coherent potential Ti{z) and the 
Green's functions G{z). Then it is easy to find that e.g. 

Ei ° = J duj 7^rr ? tij irl) Im(5 ^ (tj + i0) (8) 

andl 

Here Ei a (z) = (z — e« + \i — A^)/^f, and by definition Ga a (z) = F[T l (z)] and 

Gija(z) = jjJ2kGj:(z)e' lk ( Ri ~ Rj ) This closes the system of equations to be 
solved in a self-consistent manner. In the next section we show numerical 
calculations of the effect of correlations on the density of states in disordered 
alloy. 



3. NUMERICAL RESULTS AND DISCUSSION 

For the purpose of numerical illustration of the general approach 
sketched in previous section we calculate here the density of states (DOS) of 
perfect interacting, random noninteracting and random interacting systems. 
For the purpose of numerical analysis we shall assume U = oo limit and 
bcc crystal structure, which leads to the following canonical tight-binding 
spectrum of noninteracting carriers in clean material 

eg = — 8t cos(k x a) cos(k y a) cos(k z a). (10) 

Using t=0.0625eU leads to the noninteracting system bandwidth W = leV. 
In the following all energies and frequencies are expressed in eV. The density 
of states corresponding to this spectrum is well known. It possesses Van Hove 
singularity in the middle of the band, which extends from —8t to +8t. The 
spectrum of a clean but interacting system in the U = oo limit is replaced 
by £ 2 £g — A^ 2 ) . The corresponding density of states is thus of the same form 
but placed in the energy window [— 8i(l - n) - A (2) , +8i(l - n) - A (2) ] and 
scaled by the factor (1 — n)~ l . Note that in this U = oo limit the upper 
Hubbard band is pushed to infinity and not visible. 

The changes of the spectrum due to correlations in disordered A\_ X B X 
alloy are illustrated in figure (1) where we show the averaged D(E) and 
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Fig. 1. Comparison of the DOS of (a) disordered noninteracting Ai~ x B x 
alloy with x = 0.5, ea = 0, eb = —0.3eV with that calculated for the same 
alloy but with interacting carriers (b). Carrier density n = 0.4. 
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conditionally averaged densities of states Da(E) and Db{E) (x = 0.5, ea = 

0. £b = —0.3eV) without (Fig. 2a) and with (Fig. 2b) electron correlations. 
The carrier concentration n = 0.4. Noninteracting alloy DOS is symmetric. 
We observe strong asymmetry, the opening of the real gap in the spectrum 
of interacting carriers and the appreciable increase of the density of states 
at the fermi level (taken as E=0 in the figure). 

In conclusion, we have shown that interplay of disorder and correlations 
leads to strong renormalisation of the electron spectrum. Let us stress that 
all other properties of such systems will also be strongly affected by inter- 
action induced disorder. The calculations of dc and ac transport properties 
are in progress. 
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